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Abstract 

We consider excited lepton (1*1 or 1*1*) production in the context of effective theories with 
being four-fermion contact interaction at LHC. We also consider the two body decay mode 
of excited lepton {[*) to standard model fermion(/) and a gauge boson V{= 7, Z, 1^). 
We have performed next-to-leading order (NLO) QCD corrections to this process. In 
spite of non-renormalizable nature of the interaction, the higher order QCD corrections 
are possible and meaningful. We have shown that these corrections can be substantial 
and significant. By considering the issue of scale dependence, it is shown that the scale 
dependence of the NLO cross sections are greatly reduced as compare to leading order 
(LO) cross section. 
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1 Introduction 

The recent discovery of a Higgs like scalar particle at LHC may complete the most successful 
model in particle physics namely the Standard Model (SM). In spite of this huge success, 
there are other issues like the replication of the fermion families, dark matter, baryogenesis 
etc. that are still not understood within the framework of the SM and addressing them needs 
physics beyond the standard model. Some possible candidates are supersymmetry [1], grand 
unification [2,3] (with or without supersymmetry), family symmetries (gauged or otherwise) 
and compositeness for quarks and leptons [4] more suitable than the SM. 

The proliferation of fermion generations suggests the possibility of quark-lepton composite- 
ness. In these theories [5,6], the fundamental constituents, preons [7], experience unobserved 
strong force. At energies far above the certain scale known as composite scale A, preons are 
asymptotically free. At a composite scale A, the interaction of preons become very strong to 
form a bound state (composites) which are to be identified as quarks and leptons compositeness. 
For energies below the scale A, the effective four-fermion Lagrangian is given by [9] 



where £ and £* represents the SM and excited lepton respectively. The Lagrangian of eqn.(l) 
is not the most general the effective four-fermion Lagrangian. There can be similar operators 
involving the quarks alone or the leptons alone. However, the Lagrangian of eqn.(l) is enough 
to serve our purpose. There are other kind of four-fermion like processes mediated by a particle 
with a mass significantly higher than the energy transfer can be approximated by a contact 
interaction [9] term. The theories with extended gauge sectors, leptoquarks [11], sfermion 
exchange in a supersymmetric theory with broken i?-parity [12] etc. are the examples of such 
type of interactions. In all these cases, the heavy fields with masses Mi ^ A [13], have been 
integrated out and left with a series of such higher-dimensional terms. The terms in eqn.(l) are 
just the lowest order (in A~^) ones. However, we shall consider only terms involving qqi*i (qqii*) 
or qqi*t. 

The excited fermions can also be transformed into ordinary SM fermions through the guage 
bosons. The effective guage mediated Lagrangian [9, 10] is given by 



where G^j,, W^v and B^^ are the field strength tensor of the S'?7(3), the SU{2) and the U{1) 
gauge fields respectively. /* and / denote the excited fermion and SM fermion respectively. 
fs, f" and /' are parameters of the compositeness. Usually they are taken to be order of 1. 

It is evident that these operators may lead to significant phenomenological effects in collider 
experiments, like e'^e~ [15], eP [14] or hadronic [16,22,23]. It is quite obvious that the effects 
would be more pronounced at higher energies for a given higher-dimensional nature of C The 
best constraints on compositeness came from the Delphi [15] and CDF [16] experiments. More 
recently the measurement of the cross section [22,23] at high invariant masses set the most 
stringent limits on contact interactions of the type given in eqn.(l). 

It is well known the QCD corrections can alter the cross sections quite significantly at 
hadron colliders. As for example, a simple process like Drell-Yan, the leading order (LO) cross 
sections have been seriously underestimated. This forced us to incorporate the next-to-leading 
order (NLO) or next-to-leading log (NLL) [20,21] results in Monte Carlos codes [20] or event 
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generators such as JETRAD [18] and HERWIG [19]. Recently, contact interactions of type 
eqns.(l,2) has received much attention both CMS [22] and ATLAS [23] collaborations. They 
have searched for heavy excited lepton via channel and put the mass bound on excited 
lepton. However, there is no higher order QCD corrections existed to such heavy excited lepton 
production cross sections mediated by a contact interaction given in eqn.(l). Therefore, all 
coUider searches of contact interaction have either been based on the leading order calculations, 
or, have assumed that the higher order corrections arc exactly the same as the SM one. In this 
article, we aim to rectify this unsatisfactory state of affairs. While it may seem that the NLO 
corrections to processes driven by such non-renormalizable interactions are ill-defined, it is not 
quite true [24,27]. In particular, if the interaction can be factorized as two currents such that 
one current with colored object and other current with colored neutral object then the NLO 
QCD corrections can be done with colored current one without any difficulties. For example, 
Ref. [24] dealt with contact interaction with SM fermions. 

The rest of the article is organized as follows. In Section 2, we start by outlining the 
general methodology and follow it up with the exphcit calculation of the NLO corrections to 
the differential distribution in the dilepton (£*£, £*£*) invariant mass. Section 3 we present 
our numerical results. And finally, we summarize in Section 4. 



2 NLO corrections 



Wc consider excited lcpton(s) production in the context of contact interaction as exemplified 
by eqns.(l,2) at LHC. The processes are 
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where Pi{i = 1,2) denotes the momenta of the incoming hadrons and /j those for the outgoing 
leptons. The pj{j = 4, 5) is outgoing vector boson V{V' = 7, Z, W) momentum. Similarly, the 
momentum px carries by the inclusive hadronic state X. In the above mentioned processes, we 
have considered only two body leptonic decay of excited leptons. The excited lepton can have 
three body decay mode. This three body decay mode may have affected by QCD corrections 
which we are considering in our future work in progress [25]. The hadronic cross section is de- 
fined in terms of the partonic cross sections convoluted with the appropriate parton distribution 
functions /(f (a;) and is given by 
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where Xi is the fraction of the initial state proton's momentum carried by the parton. i.e. the 
parton momenta ki are given by ki = XiPi. For our calculational purposes, we define some 
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other variables as given below 

r = — z = — T = ZX1X2. 

S s 

As we argued above that the current-current structure of the matrix element helped the higher 
order QCD corrections possible. Therefore it is convenient to recast the matrix element for 
the abovementioned process as a sum of several current-current pieces with a "propagator" in 
between as given below 

i 

where the dots (•) denote Lorentz index contractions as appropriate and the propagators Pj are 

PvA = ^ ^^-^ 

Using this definition, the partonic cross section for the process a{ki) + b{k2) j{q) + ^{Pi) 
is given by 



where \J^"-^'^i^' denotes the square of the hadronic current, and dPS^+i the (m + l)-body 
phase space element, viz. 

dPSm-,1 = /nr((g^27r5+(pf)) (|^27r5+(g2-g^) ^^^^ 

where n = 4-|-e is the dimension of spacetime and 5~^{x) carries its usual meaning. The leptonic 
tensor, given by 

Cn'-^W ^jfl - ml)^ {2^T5^-\q -h- h)\M^^'-^^'^' \' , (11) 



is easy to calculate and which leads to 



C,,,^i.r = ( - + ^) il' = I. n (12) 



with 

c,.m = ^(0=-^-^^=^) (13) 

After substituting Cjji^ui in eqn.(9), we have the hadronic cross section, 

2^ ^(r, Q') = ^ E,y=yA PAQ') ^n'm Wff^r, Q') (14) 

where the hadronic structure function W is defined to be 

^ I, ^ ^/ " JO JO 
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The bare partonic coefficient function A is computed from the hadronic part of the matrix 
element square integrated over phase space as given below 



^{z, Q\ e) = / dPS^^, \M^^^^''\^ T,Aq): (16) 
where Tjj depends upon the spin of the current in question, viz 

TjAq) = + iJ,f-VA). (17) 

To calculate the distribution of the excited lepton pair (/*/* or /*/), one needs to compute 
the square of the hadronic matrix element p Tjj'(q), in a suitable frame so that the 

integrations over the phase space and z become easy. 

In general, the bare partonic coefficient function A{z,Q^,e) contains all three singularities 
namely ultraviolet, soft and collinear. We follow the dimensional regularization to regulate 
these singularities. The rcnormalization of the V^A-type coupling is known and can be found 
in Ref. [26,27]. We have used this rcnormalization procedure to renormalize contact interac- 
tion couplings r)ij{rilj). Defining the renormalized couplings Tiij{ri[j), we removed the ultraviolet 
singularities. By adding the contribution from the real emission diagrams [28] to the the ultra- 
violet regularized (and renormalized) virtual corrections we removed the soft singularities and 
we left with the expression containing only collinear singularities. The left over collinear sin- 
gularities is removed through the mass factorization [29]. The details of the mass factorization 
procedure can be found out [24,29]. Finally, we got the finite coefficient functions A(2;, Q^) at 
the end of the mass factorization procedure. The physical hadronic cross section is be obtained 
by convoluting these finite coefficient functions with appropriate parton distribution functions. 
By defining the renormalized parton-parton fluxes Hab{xi,X2, iJ^p) where 

Hgq{Xi,X2, Hl) = fg^(Xi,IJ,l)fp{x2, fil) + f[^(Xi,IJ,l)fg^{x2, /u|) 
H,,{X,,X2.I,1) = /f (^l,/^F)(/f (^2,/^^) + /f (X2,/^^)) (18) 

Hgg{Xl,X2,HF) = Hgg{X2,Xl,IJ,F) , 

the inclusive differential cross section can be expressed as 

= E/ / dX2 dz 5{t - ZX^X2)T''^ QvA 

dQ-^ q Jo Jo Jo 

GvA ^ Hgg{x,,x2.^J['^){^.'•;^^''\z,Q\^^^^^ 

+ {-f^</s(a;i,a^2,/^|) + i?s</(a;i,a;2,//|)}a,Ag'^^(^,//|,) 

The constants contains information of all couplings, propagators and the massive final 
state particles which is given by 

-rVA _ {ml + ml) _ {ml-mlf \ 

~ 12 AM 2Q2 2Q4 / ^^"^ 



For the vector-axial vector couplings, the results for the coefficient functions are analogous 
to the case of the SM [26] , namely 
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The SU{N) color factors in the above equations are 
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(22) 



(23) 



(24) 



3 Results and Discussion 

In the previous section, we have calculated the differential distributions with respect to the 
excited lepton (pair) invariant (either or /*/*) mass. For our interest, we have expressed 

the aforementioned differential distribution to the total cross section by integrating over Q'^. 

For the time being, we have first considered IJ^'r = l^'p = ^-iid later we have shown the scale 
dependence of our results. From the eqn.(19), the total cross section is obtained by integrating 
over and is given by 

a^^^^Ml S,A)^J (25) 

For the numerical analysis, we present our results at three different LHC energies \/S = 
7, 8, 14 TeV. Although the calculation of QCD correction does not depend on the contact 
interaction scale A, for definiteness we have used A = 6 TcV for each LHC energy unless it 
is quoted. In presenting our results, we have put all contact interaction coupling strength to 
be unity. In our numerical analysis, we have used CteqGPdf [30] parton distribution functions 
(PDFs) otherwise mentioned specifically. For Cteq6Pdf, the leading order (LO) hadronic cross 
section is obtained by convoluting the LO parton distribution (namely Cteq611) function with 
LO partonic cross section and for the NLO hadronic cross section, we have convoluted NLO 
parton distribution (namely Cteq6m) with NLO partonic cross section. We did our numerical 
analysis with Aqcd — 0.226(0.165) GeV for NLO (LO) for n/ — 5. To start with we will first 
discuss the NLO corrections of 1*1 (//*) and /*/* productions in general and later we consider 
only a particular process ll^y production where we have multiplied by the branching fraction of 
the excited lepton /* (l*) (decays to (^7)) to the production cross section 1*1 (ll*). The later 
process has been analyzed by both CMS [22] and ATLAS [23] collaborations. Although the 
two body decays of excited lepton I* (l*) (decays to ([7)) does not have any effect on QCD 
correction, for definiteness we have done this analysis so that NLO results can be read directly 
from the figures. 

In figure 1, we have plotted total cross section of excited lepton pair (1*1*) as well as one 
excited lepton and a SM lepton (/*/, II*) versus the excited lepton mass of M^, for all the light 
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M* (TeV) M, (TeV) M« (TeV) 

Figure 1: Variation of total cross-section for £*£* and £*£ {££*) production with respect to excited 
lepton mass M^, at LHC. For each set, the solid (dashed) refer to NLO (LO) cross sections. 
Upper (lower) set is represents £*£ {£*£*) for A = 6 TeV only. 

flavors ("u, (i, s-quarks). In figure 2, we have shown the variation of total cross section of one 
excited lepton and a SM lepton (/*/, //*) with respect to excited lepton mass M^, for individual 
flavor at initial state. We have done this to examine whether it is flavor dependent or not. As 
we have seen from figures 1, 2 that the cross section decreases as excited lepton mass of 
increases due to the fact that parton distribution functions fall at higher momentum fractions. 
The fall of total cross section is more lower center of mass (c.o.m.) energy \/S{= 7, 8 TeV 
say) than the higher c.o.m. energy \/S{= 14 TeV). The reason for this is higher momentum 
fraction r(= x, the Bjorken scale) and hence we are integrating over small phase space region 
at lower center of mass energy \/S. As expected, the production cross section is less and 
falls faster as compare to /*/ (//*) production cross section. All the cross sections (figures 1, 2) 
look similar behaviour except the numerical values. From figure 2, it apparently seems that the 
cross sections for different flavors at initial state are different i.e. they are flavor dependent. 
Actually it is not. This difference between the individual cross sections are due the individual 
flux difference. This implies that contact interaction is flavor blind. 

To quantify the enhancement of NLO cross section, we define a variable called X-factor 
given by 

^NLO _ _ 

Ki = ' i = total, uu, dd, ud, du (26) 

where the LO (NLO) cross sections are computed by convoluting the corresponding parton-level 
cross sections with the LO (NLO) parton distribution functions. 

In figures 3,4 we have shown the variation of iC-factor with respect to the M*. The variation 
of total i^-factor (figure 3) is about 25% — 30% for moderate values of M*(< 1 TeV) at low 
c.o.m energies (a/S* = 7, 8 TeV). At larger mass region (M* > 1 TeV), the K-factor rises very 
fast 25% — 60%. At high c.o.m energies (v^ = 14 TeV), the variation of i^-factor is about 
25% — 30% for reasonably high value of mass range (M* < 2 TeV). In figure 4, we have shown 
the K-factor variation for individual flavor only for the l*l{ll*) production process. In figures 
(3,4), the rate of fall of X-factor is much slower at higher c.o.m energy (say \/S = 14 TeV) than 
the lower c.o.m energies because of the fact that at lower c.o.m energy, we are integrating over 



6 




10^* I I I I I I I I I I I I I I I I I I I I I I I I I I I ^ 

0.4 0.8 1.2 1.6 2 0.4 0.8 1.2 1.6 2 0.4 0.8 1.2 1.6 2 

M. (TeV) M. (TeV) M. (TeV) 




0.4 0.8 1.2 1.6 2 0.4 0.8 1.2 1.6 2 0.4 0.8 1.2 1.6 2 

M„ (TeV) M, (TeV) M. (TeV) 



Figure 2: Variation of total cross-section for tI{U*) production with respect to excited lepton 
mass M^, at LHC. For each set, the solid (dashed) refer to NLO (LO) cross sections. In upper 
panel, the upper (lower) set represents uu{dd) initiated process and in lower panel, the upper 
(lower) set represents ud{du) initiated process at horn level for A = 6 TeV only. 



smaller phase space region and also higher momentum fraction. As the Bjorken x increases 
to unity, the parton distribution function falls very steeply and also the fall of LO parton 
distribution is more than the NLO parton distribution. This is the reason at lower energy, 
the K-factor increases very fast as mass M^, increases towards the center of mass energy. As 
expected (flavor independent), variation of individual K-factor is quite similar to the variation 
of total ii'-factor. One can also see from figure 4 the numerical difference between the individual 
quark K-factors is due to their respective flux difference. 



3.1 production 

The excited heavy lepton will decay into a light SM fermion and a electroweak gauge bosons 
V{= 7, Z, W) according to Eqn(2). The total NLO cross section of lepton pair (Jl) and a gauge 
boson V can be calculated by multiplying the branching ratio to Eqn.25 as given below 

a^^^^(M2,5,A) = BR{1* -> IV) J '^^^^^^^dQ\ (27) 

The partial decay widths of excited leptons for various electroweak gauge bosons are given 



by 



r(r^<v) = 3«/?^ 1-3^ 2 + ^ (28) 
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Figure 3: Variation of k-f actor with respect to excited lepton mass for A = 6 TeV at LHC. 
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Figure 4: Variation of k-f actor with respect to excited lepton mass for A = 6 TeV at LHC. 
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In figure 5, we have plotted total cross section versus invariant mass of one lepton and 
a photon M*(= Me^,Mi^)) for two PDFs CTEQ6 and MSTW 2008 [31] all the light flavors 
[u, d, s-quarks) for two different contact interaction scale A = 2, 6 TeV. The cross section 
decreases as invariant mass increases due to the fact that parton distribution functions fall 
at higher momentum fractions as mentioned above. The variation of cross section looks same 
from figure 5 for two different PDFs. In Actual practice they are not same. This can be found 
out from figures 6 and has been explained later on. From the figures 5 we see that as the 
contact interaction scale (A) increases, the cross section (both LO as well as NLO) decreases 
uniformly as A~^ as expected (from Lagrangian 1) for a fixed center of mass energy (S) and 
for different values of A, the cross section scales accordingly. Therefore, one can obtain the 
cross section (for both LO as well as NLO) for arbitrary values of A by multiplying with an 
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Figure 5: Total cross-section for iiu production at LHC. For each set, the solid (dashed) refer 
to NLO (LO) cross sections. Upper (lower) set is for A = 2(6) TeV. 

appropriate scale factor to our results. 

In figures 6, we have shown the variation of total X-factor with respect to the M*. The 
variation of iC-factor is about 25% - 35% for both PDFs CTEQ6 and MSTW 2008. The major 
difference in X-factor between the two PDFs (specially at low center of mass energy) is due 
different parameterizations of their parton distribution functions and the different data sets. In 
figures 6, the rate of fall of X-factor is much slower at higher c.o.m energy (say = 14 TeV) 
than the lower c.o.m energy as explained before. 

In our above discussion, we consider for simplest = where the cross section 

depends only on physical scales the c.o.m. energy {\/S) and the masses of final state particles 
(M*). Now we turn on another scale called factorization scale fi'p{= /i/j the renormalization 
scale). In figure 7, we have shown the factorization scale dependence. We have seen from this 
figure that the scale dependence reduced greatly at NLO cross section compare to LO cross 
section. This signifies the necessity of NLO QCD corrections. 



9 




1.18 ' ' ' ' ' ' ' ' ' ' 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Figure 6: K -factor for iiu production at three different energies of LHC. The upper (lower) set 
is for MSTW 2008 (CTEQ6) parton distribution functions. 



4 Conclusions 

To conclude, we have systematically calculated the next-to-leading order QCD corrections for 
the V ± A type contact interactions eqn.(l). Opposed to naive expectations, we have showed 
that that the QCD corrections are meaningful and reliable to such non-renormalizable theory. 

We have analyzed the variation of cross section with respect to excited lepton mass (and 
hence the invariant mass of one SM lepton and one SM gauge boson) at the LHC. The en- 
hancement of NLO cross sections over the LO cross sections expectations are found to be quite 
significant. To quantify the enhancement, we present the corresponding i^-factors in a suitable 
form for use in experimental analyzes. As it is well known that the predictions on cross sec- 
tions calculated at leading order in perturbation theory suffer scale uncertainty resulting from 
arbitrary choice of factorization scales and renormalization scales. These scale uncertainties are 
due to the absence of higher order contributions to the leading order in perturbation theory. By 
including more and more higher order contributions to theory, these scale uncertainties reduce 
gradually and the predictions become more reliable. We have showed these scale dependence 
of our results. As expected, we have seen that the scale dependences reduced greatly for the 
case of the NLO results as compared to that for the LO case. 
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